Cluster Percolation and Critical Behaviour in Spin Models and SU(N) Gauge Theories 
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The critical behaviour of several spin models can be simply described as percolation of some 
suitably defined clusters, or droplets: the onset of the geometrical transition coincides with the 
critical point and the percolation exponents are equal to the thermal exponents. It is still unknown 
whether, given a model, one can define at all the droplets. In the cases where this is possible, 
the droplet definition depends in general on the specific model at study and can be quite involved. 
We propose here a simple general definition for the droplets: they are clusters obtained by joining 
nearest-neighbour spins of the same sign with some bond probability ps, which is the minimal 
probability that still allows the existence of a percolating cluster at the critical temperature T c . By 
means of lattice Monte Carlo simulations we find that this definition indeed satisfies the conditions 
required for the droplets, for many classical spin models, discrete and continuous, both in two and 
in three dimensions. In particular, our prescription allows to describe exactly the confinement- 
deconfinement transition of SU(N) gauge theories as Polyakov loop percolation. 

PACS numbers: 64.60.Ak, 64.60.Cn, 11.15.Ha 
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I. INTRODUCTION 

The study of critical phenomena is one of the most fas- 
cinating topics in physics: phase transitions are processes 
that physicists continuously encounter in their investiga- 
tions of nature, and they can occur in an incredible vari- 
ety of systems. A phase transition is basically a change in 
the order of a system: by varying some parameter (usu- 
ally the temperature) , the microscopic constituents of the 
system, that we here call 'spins', choose a different way 
of staying together. If we go from high to low tempera- 
tures, we see that one passes from a situation in which 
the behaviour of each spin is totally independent of the 
others, to a situation in which several spins are correlated 
with each other and form ordered structures in the body 
of the system. These structures are local realizations of 
the new phase Q|. If we further lower the temperature, 
the size of the ordered regions increases, until most of 
them stick to each other and form a domain which spans 
the whole system. When this happens, the system is in 
a new state of order, i.e. in a new phase. 

This mechanism led already at the end of the 40 's 
to the idea that the physics of the phase transition is 
basically governed by the ordered domains built by the 
spin-spin correlations, and not by the individual spins. 
If the degrees of freedom relevant for the phase change 
are the ones of sets of particles, and not of single parti- 
cles, it is likely that they do not depend on the details of 
the microscopic interaction, but only on its gross features 
(e.g. symmetries). That could explain why whole classes 
of systems, ruled by dynamics which appear very differ- 
ent from each other, happen to have the same behaviour 
at the phase transition. In particular, for second-order 
phase transitions, this simple picture could allow to ex- 
plain the universality of the critical indices. 

Percolation theory || is the ideal framework for a ge- 
ometrical description of phase transitions. The percola- 
tion phenomenon takes place when geometrical clusters, 
formed by elementary objects of some system, stick to 



each other giving rise to an infinite network, that spans 
the whole system. The analogy with the phase transi- 
tion mechanism described above is evident, but there is 
much more than that. The real amazing thing is the fact 
that, in spite of the apparently different nature of per- 
colation, which is a purely geometrical phenomenon, and 
second-order thermal phase transitions, these two types 
of processes are formally identical, as they are character- 
ized by the same basic features, which are: 

• Power law behaviour of the variables near the crit- 
ical point, with relative exponents; 

• Scaling relations between the exponents; 

• Universality of the critical indices. 

The temptation to identify a continuous thermal phase 
transition with a simple percolation transition is then 
very strong. This identification is possible provided one 
can establish a one-to-one correspondence between ther- 
mal and percolation variables. The main percolation 
variables are: 

• the percolation strength P, i.e. the probability that 
a site chosen at random belongs to a percolating 
cluster; 

• the average cluster size S, 



(1) 



where n s is the number of clusters with s sites and 
the sums exclude eventual percolating clusters. 

Suppose we have defined how to group the spins of the 
configurations of a given model in clusters. Such clusters 
are the 'physical' clusters or droplets of the model if the 
following conditions are satisfied: 
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• the percolation point coincides with the thermal 
critical point; 

• the connectedness length (average cluster radius) 
diverges as the thermal correlation length (same 
exponent); 

• the percolation strength P near the threshold varies 
like the order parameter m of the model (same ex- 
ponent); 

• the average cluster size S diverges as the physical 
susceptibility x (same exponent). 

The hrst studies in this direction concerned the Ising 
model. The simplest clusters one can think of are just 
the magnetic domains, i.e. the clusters built by join- 
ing nearest-neighbour spins of the same sign. In two di- 
mensions these clusters indeed percolate at the critical 
temperature of the Ising model [Q; however, the critical 
percolation exponents differ from the thermal ones ||. 
So, the magnetic domains are not the Ising droplets. It 
was soon realized that the problem of these clusters is 
that they are too big due to purely geometrical effects. 
The average cluster size would be non-negligible also at 
extremely high temperatures, just because spins of the 
same sign may happen to lie close to each other, even 
in the absence of a physical correlation between them. 
In order to eliminate these artificial geometrical correla- 
tions, one can introduce a bond probability ps and join 
nearest-neighbouring spins with this probability, which 
automatically reduces the size of the clusters. Since the 
correlation changes with the temperature T, ps must be 
as well a function of T. If one chooses the expression 
Pb = 1 — exp(— 2J/kT) || (J is the Ising coupling), the 
corresponding site-bond clusters are indeed the critical 
droplets of the Ising model, in any dimension M |) . 

Nobody can yet say whether the phase transition of 
every model can be geometrically described as a perco- 
lation transition, i.e. whether one can always define the 
droplets. The Ising result, which is valid more in gen- 
eral for the g-state Potts model, can easily be extended 
to several spin systems, both discrete and continuous 
§, 0, 0. In general, one finds that each interaction 
between a pair of spins corresponds to a bond in the per- 
colation picture with an analogous bond probability as 
in Ising. In models with several spin-spin interactions of 
the same type (all ferro- or antiferromagnetic), one can 
still define a percolation picture Q by putting bonds 
between any pair of interacting spins with some prob- 
ability, but the picture becomes quite involved: if two 
interacting spins are far from each other, the geometrical 
bond between them looks virtual, as the two spins are 
geometrically disconnected. 

Recent investigations [|l2[ aimed at recovering the 
importance of the role of geometrical connectivity in the 
mapping between percolation and thermal critical be- 
haviour. In [jl3| it was shown that, for a wide class of 
bidimensional models, one can define simple site-bond 
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clusters which show all features the droplets should have. 
For the models where a rigorous mapping between per- 
colation and critical behaviour is possible, such site-bond 
clusters are in general different from the "exact" droplets 
(see |,fuj 0), which are in general more complex, even 
if their behaviour at criticality is identical. Moreover, the 
result remains valid as well for models with competitive 
interactions (e.g. ferromagnetic + antiferromagnetic), for 
which an exact definition of the droplets is, at present, 
missing. In fl3|| one examined theories with center sym- 
metry Z(2) and Z(3), such that their critical behaviour 
is in the universality class of the model obtained by re- 
moving all interactions except the nearest-neighbour one 
(Ising for Z(2), 3-state Potts for Z(3)). Therefore, for 
these models the nearest-neighbour spin-spin coupling is 
the fundamental interaction which determines the be- 
haviour at the phase transition. This is probably the 
reason why, if one just considers geometrical connections 
between nearest-neighbours, weighted by some suitable 
bond probability, the corresponding clusters are at least a 
good approximation of the critical droplets of the model. 
The previous argument is of course independent of the 
number d of space dimensions of the system. For this 
reason we believe that the result of [ 
general, i.e. for d > 2 as well. In t 
the 2-dimensional results presented in |13| and provide 
numerical evidence, based on Monte Carlo simulations, 
that our droplet definition holds true also in three dimen- 
sions. 

We stress that the original target of our investiga- 
tions was to provide a geometrical description in terms 
of percolation of the confinement-deconfinement transi- 
tion in SU(N) gauge theories. Early attempts focused 
on SU(2) pure gauge theory, whose deconfining transi- 
tion is second-order and in the universality class of the 
Ising model Q . The strategy we followed at that stage 
was to approximate the gauge model by means of sim- 
pler Polyakov loop effective theories for which an ex- 
act droplet definition exists 15 [l(|. In this way one 
finds only an approximate solution of the problem, which 
strongly depends on the specific lattice regularization 
one chooses |l7j]. This seemed to us unsatisfactory: the 
droplet prescription we propose here solves the problem 
in a simple and general way. 

The paper is divided as follows: in Section II we dis- 
cuss more in detail our droplet definition; in Section III 
we present the results of our simulations, distinguishing 
between spin models and SU (2) pure gauge theory; fi- 
nally the conclusions of our work are exposed. 



II. THE DROPLET DEFINITION 

Our droplet candidates are clusters built by joining 
nearest-neighbour spins of the same sign with a bond 
probability ps- they are then uniquely defined once we 
specify ps- In Jl3| , the following criterion was proposed: 
is the minimal probability that still makes percolation 
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possible at the critical temperature T c . This special min- 
imal probability, that we will call pck , depends on the 
model at study. We recall that pck is in general a func- 
tion of the temperature, like in the Ising model (where 
Pck{T) = 1 — exp(— 2J/kT)). We are only interested 
in the behaviour at the transition, i.e. near T c . There- 
fore what matters is basically only the value pck(T c ) of 
the bond weight at T c ]is| . This is why, referring to pp 
(Pck), we use the term "value" instead of " expression" . 

In two dimensions, for any ps < Pck, the percola- 
tion temperature T p < T c and the exponents are in the 
universality class of 2D pure random percolation. On 
the other hand, for any pg > pck (including the pure- 
site case pb — 1), T p — T c , but the exponents do not 
coincide either with the thermal or with the random per- 
colation ones. For pg = pck and only in this case, the 
site-bond clusters satisfy all conditions required for the 
droplets, i.e. both the critical temperature and the ex- 
ponents of the geometrical transition are equal to the 
thermal counterparts. We see that there is a whole range 
of pb values, i.e. pck < Pb < 1, for which the relative 
site-bond clusters begin to percolate exactly at the onset 
of the thermal transition. Such feature is specific of bidi- 
mensional lattices. In three dimensions we shall see that 
there is just a single value of the bond probability pb for 
which the percolation threshold is exactly at T c . If pb 
is greater than this value, the clusters begin to percolate 
at some T p > T c . However, this asymmetry between the 
2D and the 3D cases does not represent a serious prob- 
lem. In the Ising model, by using the Fortuin-Kasteleyn 
bond weight pck = 1 — exp(— 2J/kT), one obtains the 
correct droplets in any dimension |7], M. In particular, 
for the 3-dimensional Ising model, pck{T c ) is necessarily 
the unique ps value for which the two thresholds coincide 
p9| . This special ps value is again the minimal proba- 
bility one needs in order to have a percolating cluster 
at T c , so that the "criterion of the minimal bond proba- 
bility" adopted in |l3| would lead to the correct droplet 
definition in three dimensions too, at least in the Ising 
case. We then assume that such criterion is valid more 
in general, and we shall verify its validity by computer 
simulations of two models in three dimensions, the 0{2) 
spin model and SU{2) lattice gauge theory. 



III. RESULTS 
A. Numerical Analysis 



is determined, all lattice sites are grouped in clusters by 
means of the algorithm devised by Hoshen and Kopel- 
man p0| ; for the cluster labeling we have always used 
free boundary conditions. After that we are left with a 
set of clusters of various sizes, and we can calculate the 
percolation variables. If a cluster connects the top with 
the bottom side (face in 3D) of the lattice, we say that it 
percolates |2l[. Besides the percolation strength P and 
the average cluster size S, we also calculate the size Sm 
of the largest cluster of the configuration, since from it 
one can derive the fractal dimension D of the spanning 
cluster at the threshold |22j . At each iteration the energy 
density e and the lattice average m of the order parame- 
ter of the thermal transition were also stored Q . 

The first step is of course the determination of the 
percolation temperature. This can be effectively done 
by using a variable that can be extracted from the data 
sample of the percolation strength P. Suppose we have 
performed a number Nj of iterations for one of our mod- 
els at a given temperature and lattice size. Looking at 
the column of the P data, say Np the number of configu- 
rations of our sample with (at least) a percolating cluster 
(for those configurations P^O). The ratio Np/N] is the 
percolation cumulant II, which shares the same properties 
of the well-known Binder cumulant, namely: 

1. if one plots II as a function of T, all curves corre- 
sponding to different lattice sizes cross at the same 
temperature T p , which marks the threshold of the 
percolation transition; 

2. the percolation cumulants for different values of the 
lattice size L coincide, if considered as functions of 
tpL 1 /^ (t p = (T — T p )/T p , v v is the exponent of 
the connectedness length); 

3. the value of II at T p is a universal quantity, i.e. it 
labels a well defined set of critical indices. 

The first property suggests that it is enough to make 
simulations on two different lattices to determine the crit- 
ical point. The result is of course the more precise the 
larger the size of the lattices. 

After evaluating the percolation temperature T p , we 
extracted the values of the critical indices by means 
of standard finite-size scaling techniques at the critical 
point. If corrections to scaling do not play an impor- 
tant role, the finite-size scaling laws of the percolation 
variables at T p are given by the simple formulas 



Our aim is to investigate the percolation transition of 
special site-bond clusters, determining in particular the 
percolation temperature and the critical exponents. 

To produce the equilibrium configurations we made use 
of standard Monte Carlo algorithms, like Metropolis or 
heat bath; for some models we adopted cluster updates, 
like the Wolff algorithm, which allows to reduce sensibly 
the correlation of the data and save a lot of CPU time. 
At each iteration, once the configuration to be examined 



P{T p ) oc Z,-V"»> (2) 
S{T p ) oc Ui*/"> (3) 
S M (T P ) cx L D , (4) 

where L is the lattice side and /3 P , 7 P are the expo- 
nents that rule the power law behaviour at criticality of 
P and S, respectively. In order to improve the precision 
of the fits and to keep disturbing finite-size effects under 
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control, for each model four to six different lattices were 
used. In all our analyses we found that corrections to 
scaling do not perturb appreciably the leading behaviour 
expressed by Eqs. (2)-(4). 

B. Spin Models 

We start by reviewing the investigations of the bidi- 
mensional models presented in ||l3l| . There, we analyzed 
two classes of systems: models with Z(2) global symme- 
try and a magnetization transition with Ising exponents 
and models with Z(3) global symmetry and a magne- 
tization transition with exponents belonging to the 2- 
dimensional 3-state Potts model universality class. The 
spin systems belonging to the first group are: 

1. the Ising model, TL = — JJ2ij s i s j (J > 0, Sj = 

±i); 

2. a model with nearest-neighbour (NN) ferromag- 
netic coupling and a weaker next-to-nearest (NTN) 
antifcrromagnetic coupling: TL = — Jj J2nn s i s 3 ~ 
J2T,NTN s i s J ( J i > 0, J 2 < 0, \J 2 /Ji\ = 1/10, 

Si = ±1); 

3. the continuous Ising model, TL = —J^ZijSiSj 
(J>0, -l<Si< + l). 

The models belonging to the second group are: 

1. the 3-state Potts model, TL = —JJ2ijH s ii s j) 
(J>0, Si = 1,2,3); 

2. a model obtained by adding to 1) a weaker next- 
to-nearest (NTN) antiferromagnetic coupling: TL = 

-JiY,NN S ( s i> s i) ~ j zI2ntn S ( s i' s j) ( J i > °i 
J 2 < 0, \ J 2 /Ji\ = 1/10, Si = 1,2,3); 

The strategy we followed in the numerical analysis was 
to tune by hand the value of the bond probability ps un- 
til, at the critical temperature T c of the model, the per- 
colation cumulant II takes the same value for each of the 
lattices we used. The smallest pb value for which this is 
still possible is the minimal bond probability pck we look 
for. The threshold value of II gives a strong indication on 
the universality class of the geometrical transition (prop- 
erty 3 of the percolation cumulant). To calculate the 
error on pck we decreased ps until the II values of all 
lattices at T c were offset by more than one a. 

From a strictly numerical point of view one should take 
care to interpret the data of the simulations when ps 
is close to pck- In this case, in fact, the system finds 
itself in the neighbourhood of a discontinuity and if the 
lattice is not large enough, its behaviour is influenced by 
that. The simulations on small lattices would produce 
configurations which represent a sort of mixture of the 
two situations at pb = Pck and pb^Pck- To recover 
the real behaviour of the system one should then go to 
very large lattices and disregard the small ones. 



In two dimensions, as we said above, there is a whole 
range of ps values such that the onset of the percolation 
transition is exactly at the thermal critical point. The 
calculation of the exponents is then crucial to distinguish 
the various geometrical transitions. 

We start by discussing the spin systems with Z(2) sym- 
metry. Fig. [j] shows the threshold value of the percola- 
tion cumulant II in the Ising model as a function of the 
bond probability ps- For ps ~ Pck we plotted the value 
of II for the largest lattice we took (1000 2 ), as IT changes 
sensibly with the lattice size for the reason we explained 
above |24j. 





2D Z(2) SP Univ. Class- 




2D Ising Univ. Class 




2D RP Univ. Class 







0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



Bond Probability p B 

FIG. 1: 2D Ising model: variation of II at the percolation 
temperature T p with the bond weight pb- 

The vertical dashed line in the plot marks the minimal 
probability p C K = 1 - exp(-2 J/kT c ) = 0.58578. We see 
that the cumulant is quite stable to the right and to the 
left of pck , and that the plateau values correspond to 
two different universality classes. For pb < Pck the site- 
bond clusters percolate at a temperature T p < T c , and 
the percolation exponents are in the random percolation 
universality class (25|. When ps > Pck T p — T c and the 
exponents belong to a special universality class (we call 
it 2D Z(2) site percolation universality class because it 
includes the pure site percolation case, Pb = 1, and is 
the same for all Z{2) models we studied). For »b = Pck 
we recover the Fortuin-Kasteleyn mapping Q and the 
site-bond clusters are the exact critical droplets of the 
system. We see that in this case the II threshold value 
does not lie on either of the plateaus, since the exponents 
are now in a different universality class, i.e. in the class 
of the 2D Ising model. 

For the other Z(2) spin systems the situation is anal- 
ogous and we could show identical pictures as Fig. [l|, 
except that the minimal probabilities pck arc different 
(p CK = 0.583(1) for Model 2 and p CK = 0.6115(9) for 
Model 3). We remark that for the continuous Ising model 
we bound nearest-neighbour spins of the same sign, in- 
dependently of their absolute values, although they play 
a key role in the definition of the "exact" droplets M. 
However, there is no proof of the existence of an exact 
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correspondence between the percolation transition of the 
site-bond clusters for pb = Pck and the magnetization 
transition. So it is essential to calculate precisely the crit- 
ical exponents to show that the "minimal" clusters are 
indeed droplets for the system. The results are shown in 
Table |, where we can see that the agreement with the 
critical indices of the 2D Ising droplets is very good. 









D 
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2D Ising 


1/8=0.125 


7/4=1.75 


15/8=1.875 


0.585(1) 


Model 2 


0.131(10) 


1.742(12) 


1.862(20) 


0.583(4) 


Model 3 


0.121(9) 


1.764(14) 


1.870(11) 


0.587(3) 



TABLE I: Critical percolation indices for the site-bond clus- 
ters of the Z(2) models when pb = Pck, compared with the 
values of the 2D Ising droplets. 
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FIG. 2: 2D 3-state Potts model: variation of If at the perco- 
lation temperature T p with the bond weight pb- 
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0.932(2) 



TABLE II: Critical indices of the percolation transition of 
site-bond clusters when ps > Pck, for the two groups of spin 
systems we considered. 



prove that the minimal clusters are droplets for the sys- 
tem we can only rely on the numerical evaluation of the 
exponents. We report our estimates in Table 111: the 
agreement with the indices of the 3-state Potts droplets 
is good. 
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0.649(9) 
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0.143(17) 


1.725(21) 


1.858(18) 


0.646(11) 



TABLE III: Critical percolation indices for the site-bond clus- 
ters of the Z(3) models when ps = Pck- 



In three dimensions, as we said in Section II, the sit- 
uation looks quite different. We performed simulations 
of the 3D Ising model, using several values for the bond 
probability above and below pck = 1 — exp(— 2 J/kT c ) = 
0.35808... We studied the variation of II with ps as we 
did in two dimensions: the result is shown in Fig. |^. For 
Pb Pck, the percolation temperature T p is different 
from the magnetization temperature T c and we recover 
the 3D random percolation exponents (see ^5|). Only 
for pb = Pck is T p = T c and we could eventually get 
the thermal critical indices (for the 3D Ising model this 
is exactly what happens). 



As far as the analysis of the Z(3) spin systems is con- 
cerned, we attain the same conclusions. Fig. || shows 
the dependence on ps of II at the percolation threshold 
for the 3-state Potts model. We notice that we obtain 
the same pattern we found for the Z(2) models (see Fig. 
|l|) . Since the magnetization transition is characterized by 
exponents which are different from the Ising ones, one ex- 
pects to find another set of critical indices to the right of 
Pck = 0.6339736.... As we can see from the figure, the 
plateau for ps > Pck lies indeed at a different height 
compared to the one we have in the Ising plot (0.932 in- 
stead of 0.9832). Analogously as we did above, we call 
this new set of critical indices 2D Z(3) site percolation 
universality class. The values of these indices and of their 
Z(2) counterparts were predicted in f2fj| , ^7j| : our numer- 
ical findings confirm such theoretical predictions, which 
are listed in Table || 

For the Z(3) model with competitive interactions the 
results are the same, apart from the value of the min- 
imal probability (pck = 0.61(1)). Again, in order to 
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FIG. 3: 3D Ising model: variation of II at the percolation 
temperature T p with the bond weight ps- 

We notice that, even if the 3D pattern is very different 
from the 2D one, the special (Fortuin-Kasteleyn) proba- 
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bility pck is still the smallest probability for which the 
site-bond clusters can form a percolating cluster at T c 
(for pb < Pck, T p < T c and all clusters at T c are finite). 

We want to verify whether this is also valid for other 
3D spin systems, and we analyze here the 0(2), or XY, 
model. For 0(n) models a rigorous mapping between 
percolation and critical behaviour was established in [jll| . 
The droplets are built in two steps: 

1. choose a random vector r of 0(n); 

2. bind together any pair of nearest-neighbouring 
spins Sj, sj with the probability 
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p{i,j) = 1 — exp{mm[0, — 2/3(si-r)(sj-r)]} (5) 
09 = J/kT). 

Such droplets are just the clusters devised by Wolff 
in his famous algorithm j28| for 0(n) spin systems. We 
see that only pairs of spin vectors having both a posi- 
tive/negative projection on the random vector r can be 
joined to each other. The random vector r, therefore, di- 
vides the spin space in two hemispheres, separating the 
spins which have a positive projection onto it from the 
ones which have a negative projection. The droplets are 
made out of spins which all lie either in the one or in the 
other hemisphere. In this respect, we can again speak 
of 'up' and 'down' spins, like for the Ising model. In 
addition to that, the bond probability is local, since it 
explicitly depends on the spin vectors Si and sj, and not 
only on the temperature like the Fortuin-Kasteleyn fac- 
tor. 

The situation is similar as in the 2D continuous Ising 
model we considered above, and we proceeded in the 
same way, i.e. we reduced the 0(2) configurations to 
Ising configurations, according to the sign of the projec- 
tion of the spins on r, so disregarding the length of the 
projection. The bond weight we introduced is the 
same for each pair of nearest-neighbouring sites. 

For our simulations we applied the Wolff algorithm and 
used four lattices: 24 3 , 48 3 , 72 3 and 96 3 . At each run, 
40000 to 100000 measurements were taken. We tuned 
the bond probability so to make the percolation point 
coincide with the magnetization point. At the end we 
found the same scenario that we had seen for the Ising 
model (Fig. §). 

Our estimate of the minimal bond probability is pck — 
0.374(1). Finally we calculated the critical indices of the 
percolation transition when pb = Pck- Fig. [| shows the 
corresponding finite size scaling plots of the percolation 
strenght P (top) and the average cluster size S (bottom) 
at the critical point. The \ 2 of the fits improves consid- 
erably if the smallest lattice (24 3 ) is excluded, this is why 
we put only three points in the plots. 

The critical indices we extracted agree with the ther- 
mal 0(2) values (Table |fv|). 



FIG. 4: 3D 0(2) model: variation of II at the percolation 
temperature T v with the bond weight pb- 



100 



Lattice Side L 



10000 
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FIG. 5: 3D 0(2) model: log-log plots of P (top) and 5" 
(bottom) at T c versus the lattice side L for the minimal site- 
bond clusters. 



C. SU(2) Pure Gauge Theory 

The deconfining transition from hadronic matter to a 
plasma of quarks and gluons has been object of inten- 
sive investigations over the last two decades. Though 
the concrete goal is to try to produce the quark-gluon 
plasma by means of high energy heavy ion collisions in 
the lab, i.e. in tiny and in general non-equilibrated fire- 
balls, it is crucial from a theoretical point of view to study 
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D 


3D 0(2) 


0.5189(3) 


1.9619(5) 


2.4808(8) 


Perc. Exponents 


0.530(15) 


1.971(13) 


2.484(7) 



TABLE IV: Critical percolation indices for the site-bond clus- 
ters of 3D 0(2) when pb = Pck\ for c omp arison we also 
report the 0(2) thermal exponents (from |29[ ) . 



with 



L = \(L S ) 



(8) 



1 Nt 



(9) 



the ideal situation of an infinite system of strongly inter- 
acting matter in thermal equilibrium at a temperature 
T. This could be effectively done after the discovery of 
the lattice approach [p0| , and indeed finite temperature 
Quantum Chromodynamics (QCD) has been extensively 
simulated on the lattice since then. 

The group that rules the gauge invariance of QCD is 
SU(3), which is non-abelian. Because of that the gauge 
fields are self-interacting and it makes sense to study sys- 
tems constitued only by gluons. This simpler situation 
is described by the so-called SU(3) pure gauge theory. 
Since any SU (N) group is non-abelian, the study of the 
relative pure gauge theories may be of interest also for 

Suppose we have a d-dimensional box containing glu- 
ons at a temperature T. The discretization of space-time 
returns a (d + l)-dimensional lattice, with N a spacings 
in each space direction and N T spacings in the imaginary 
time (or temperature) direction. The partition function 
of finite temperature SU(N) pure gauge theories on this 
lattice takes the form 



Z(N a ,N T ;g 2 ) = I [J dUy exp[-S(U)}, (6) 

links 

where S(U) is the Wilson action 

s ( u ) = ™ E i 1 - ^ Re Tr uuuu ) • ( ? ) 

Here g is the (temperature-dependent) coupling and 
Uij the so-called link variable, which is a function of the 
gauge fields set between a pair of nearest-neighbouring 
sites i and j. The product in Eq. runs over all links 
of the lattice, the sum in Eq. (Q) over all the smallest 
closed paths (plaquettes), which are formed by four links; 
UUUU is the product of the link variables corresponding 
to each side of a plaquette. 

All SU(N) pure gauge theories undergo a transition 
from a phase in which the gluons are bound in glueballs to 
a phase of free gluons. Such deconfining transition is due 
to the spontaneous breaking of a global Z(N) symmetry 
which results from the periodicity of the gauge fields in 
the temperature direction ]3l[ |. The order parameter is 
the lattice average of the Polyakov loop, defined as 



The product in (0) runs over all the U's in the tem- 
perature direction taken at a given spatial site x. In the 
confined phase L = 0, whereas at deconfinement L =/= 0. 
The main features of the deconfining transition are then 
all in the Polyakov loop configurations one obtains by 
projecting out the temperature direction of the lattice 
through the matrix product of Eq. (||). 

There are conjectures suggesting that the deconfining 
transition of SU(N) pure gauge theories is intimately 
related to the magnetization transition of the iV-state 
Potts model, with which they share the Z(N) symmetry 
p2| . In particular, it was predicted that if the deconfin- 
ing transition is second-order, the critical indices are in 
the universality class of the corresponding Potts model: 
this prediction has been confirmed by computer simu- 
lations without exceptions. This is actually the reason 
why we investigated simple Potts-like spin systems. As 
our droplet prescription seems to work for these models, 
we tried to see whether it is correct for SU (N) gauge the- 
ories as well, at least in the cases in which the deconfining 
transition is continuous. In this way we would have for 
the confinement-deconfinement transition the same geo- 
metrical picture as for magnetization in the Potts model. 
If we take a typical Polyakov loop configuration of an 
SU(N) theory at a certain temperature, there will be 
areas where L takes negative values, and areas where L 
takes positive values. Both the positive and the negative 
"islands" can be seen as local regions of deconfinement. 
As long as there are finite islands of both signs, decon- 
finement remains a local phenomenon and the whole sys- 
tem is in the confined phase. When one of this islands 
percolates, i.e. it becomes infinite, we can talk of decon- 
finement as a global phase of the system. 

We examined the simple case of SU(2) pure gauge the- 
ory, in two and three space dimensions. In both cases we 
focused on one and the same lattice regularization, cor- 
responding to N T = 2; however, the result should hold 
for any N T (see Conclusions). We stress that for SU(2) 
the Polyakov loop L% is a real number, which can take 
all possible values in a range (from -1 to +1 according to 
our normalization choice). We are in the same situation 
as in the continuous Ising model, and again we will con- 
sider only the sign of Lg, disregarding its absolute value. 
So, the Polyakov loop configurations reduce themselves 
to Ising-like configurations, that we treated in the same 
way as in the previous section. 

We start by exposing the results for 2 + 1 SU(2). 
We carried on our simulations on four different lattices, 
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64 2 x2, 96 2 x2, 128 2 x2 and 200 2 x2. For each run we col- 
lected from 10000 to 40000 measurements. The result of 
our analysis is identical as for the Z(2) spin systems (Fig. 
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FIG. 6: 2 + 1 SU(2): variation of II at the percolation tem- 
perature T p with the bond weight pb- 



If we compare Fig. || to Fig. |] we see no differences, 
except in the value of the minimal bond probability (for 
SU(2), pck — 0.6275(5)). 

This is truly remarkable and shows the substantial 
analogy of the two cases. The exponents we calculated 
for pb = Pck are in good accord with the 2D Ising 
ones, which shows that the minimal clusters are critical 
droplets for SU(2) as well (Table 0). 





Pp/Vp 


Jp/Vp 


D 


n at T p 


2D Ising 


1/8=0.125 


7/4=1.75 


15/8=1.875 


0.585(1) 


SU(2) Perc. 


0.140(19) 


1.761(17) 


1.882(18) 


0.586(5) 



TABLE V: Critical percolation indices for the site-bond clus- 
ters of 2+ 1 SU{2) (N T = 2) when ps = Pck, compared with 
the values of the 2D Ising droplets. 



Let us now turn to 3 + 1 SU(2). We used again four 
lattices, namely 12 3 x2, 20 3 x2, 30 3 x2 and 40 3 x2. The 
number of measurements ranged from 20000 to 50000 for 
each run. The pattern we derived for LI at the percolation 
threshold (Fig. [j]) is the same as the one of the 3D Ising 
model (see Fig. ||), only the value of pck is different 
(here is pck = 0.385(4)). 

We notice in particular that for ps = Pck the thresh- 
old value of LI sits at the height corresponding to the 
3D Ising universality class. The exponents relative to 
the transition of the minimal clusters agree within errors 
with the 3D Ising exponents, as we expected (Table VI). 
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FIG. 7: 3 + 1 5(7(2): variation of II at the percolation tem- 
perature Tp with the bond weight ps- 





Pp/Vp 


Jp/Vp 


D 


n at T p 


3D Ising 


0.5187(14) 


1.963(7) 


2.4817(21) 


0.423(5) 


SU(2) Perc. 


0.54(3) 


1.962(8) 


2.477(9) 


0.424(5) 



TABLE VI: Critical percolation indices for the site-bond clus- 
ters of 3+ 1 SU{2) (N T = 2) when ps = Pclk, compared with 
the values of the 3D Ising droplets (from ]33|). 



IV. CONCLUSIONS 

We have shown that the phase transition of many 
systems can be interpreted as a geometrical percolation 
transition of simple site-bond clusters, like in the Ising 
model. We have found a general criterion to define the 
"correct" bond probability pck of geometrical connec- 
tion between nearest-neighbouring sites carrying spins of 
the same sign: pck is the minimal probability for which 
it is still possible to have a percolating cluster at the 
critical temperature T c of the system. We remark that 
this criterion looks a bit artificial, since it imposes by 
hand the coincidence of the percolation with the thermal 
threshold. So, if one studies a new model, it would be im- 
possible to define the droplets until one finds the critical 
temperature of the system. However, this is not relevant 
for us, as our aim was just to show that the droplets ex- 
ist. We investigated a wide variety of models, from spin 
systems (discrete and continuous) to SU(2) lattice gauge 
theory, both in two and in three dimensions. To the ex- 
tent of these models our recipe provides indeed a solution 
of the problem, valid even for systems for which "ex- 
act" droplets could not so far be identified (models with 
competitive interactions, SU(2)). The generality of the 
solution is clearly shown by its validity for a complex the- 
ory like SU(2): the effective theory of the Polyakov loop 
for SU (JV) gauge theory consists of a mixture of many 
different interactions, short- and long-ranged, ferromag- 
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netic and antiferromagnetic, including couplings between 
more than two spins (like plaquette-interactions, six-spin 
couplings, etc.) and self-interactions. Our analysis is en- 
tirely numerical, but the result is most likely exact, as it 
is for Ising. 

The droplet definition we propose puts in evidence the 
key role of geometrical connectivity in the mechanism of 
the phase transition. This is easy to understand for the 
models we have considered, where the nearest-neighbour 
coupling is by far the most important compared to even- 
tual longer-ranged interactions and determines the criti- 
cal behaviour. However this may not be valid when the 
strength of other couplings is comparable to the nearest- 
neighbour one: in this case the phase transition could 
be influenced as well by the other interactions (e.g. the 
critical indices might change) and a droplet definition 
based only on nearest-neighbour connections is probably 
inadequate. 

In our analysis of the models with Z(2) symmetry we 
have seen that the essential spin feature for the droplet 
definition is the Z(2) variable, i.e. the sign of the spin. 
In fact, we built the clusters in the same way, no mat- 
ter if the spin is discrete like in Ising or continuous like 
in SU(2). This shows once more the crucial role played 
by the Z(2) symmetry, whose spontaneous breaking is 
indeed responsible for the phase transition. Since the re- 
lationship between SU(N) pure gauge theories and Z(N) 
spin models holds for any N, our result should be valid 
for all SU (N) pure gauge theories with a continuous de- 
confining transition, i.e. also for SU(3) and SU(4) in 
2 + 1 dimensions. The special lattice regularization of 
the gauge theory does not play a role, as the critical be- 
haviour has the same features for any N T . So, we found 
that the confincmcnt-deconfinement phase transition of 
SU(N) pure gauge theory, if second-order, is equivalent 
to a percolation transition of special site-bond clusters 
of like-signed Polyakov loops: this is the most important 



result of our work. Very recent studies show that perco- 
lation can help to describe as well the chiral transition of 
special fermion lattice models [[34] . 

We focused only on systems with a continuous second- 
order transition because the percolation phenomenon is 
typically smooth. This does not mean that we cannot use 
percolation to describe discontinuous phase changes. The 
Fortuin-Kasteleyn mapping Q is valid for any q— state 
Potts model in any dimension, and the equivalence be- 
tween the magnetization and the percolation order pa- 
rameter holds for first-order phase transitions as well. In 
this case at the critical point there is a coexistence of 
a paramagnetic and a ferromagnetic phase, which corre- 
spond to two different "geometrical phases", character- 
ized by small and large clusters, respectively. Therefore, 
the percolation order parameter jumps at the threshold 
exactly as the magnetization (from zero to a non-zero 
value). Because of that we believe that the criterion of 
the minimal bond probability we adopted here can be 
directly extended to systems undergoing a discontinu- 
ous phase transition. Now it is however more difficult to 
support the idea of droplets because there are no criti- 
cal indices to reproduce. One could nevertheless try to 
see whether the percolation variables (P, S) vary near 
T c like their thermal counterparts (m, x)- 111 this way 
one could study in particular the "physical" confincmcnt- 
deconfinement transition of SU (3) pure gauge theory in 
3+1 dimensions (quenched QCD). 
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